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In  [1],  we  gave  another  proof  of  the  boundedness  of  L2~projections  onto 

splines  on  a  geometric  mesh.  In  this  paper,  we  obtain  the  sharp  lower  bound 
for  the  inverse  of  the  corresponding  B-spline  Gram  matrix.  I.e. 


I 


00 


n2k-1(qr>q) 

n2k-1(_qr,<I) 


>  2k- 1 , 


for  r*k,  k-1. 
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SIGNIFICANCE  AND  EXPLANATION 


-  Least-squares  approximation  by  polynomial  aplinea  is  a  vary  effective 
means  of  approximation/  particularly  when  the  knots  are  appropriately 
nonuniformly  spaced  to  adapt  to  the  particular  behaviour  of  the  function  being 
approximated.  Unfortunately,  the  stability  of  this  process  has  been 
established  only  for  nearly  uniform  knot  sequences.  The  stability  can  be 
linked  to  the  norm  of  the  inverse  of  the  Gram  matrix  of  a  (appropriately 
scaled)  B-spline  basis.  In  an  earlier  report  ]£Y,  we  studied  an  important 
special  case,  that  of  a  geometric  knot  sequence  and  there  showed  the  norm  of 
the  inverse  of  that  Gramian  to  be  bounded  independent  of  the  mesh  ratio.  -s 
In  the  present  report,  we  continue  these  investigations  and  show,  in 
particular,  the  surprising  fact  that  the  norm  of  the  inverse  of  the  Gramian  is 
least  (i.e.,  the  stability  is  greatest)  when  the  mesh  is  most  nonuniform. 


/ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


L^-Lower  Bound  of  Lj -Pro j actions  onto  Splines 


on  a  Geometric  Meah 


Y.  Y.  Feng*  and  J.  Xozak** 


1 .  Introduction 


He  begin  with  the  explanation  of  some  notation*. 

j  n  .  n  j 

K  (Xiqi  t  *  -  J  <-ln  (  )  U  (q^-X)  9  the  generalized  Euler-Frobenlus  polynomial  of 

i  j-0 

jfi 


nit  i-0 


order  n.  t  i  »  in  q 

(!) - n‘ 


>  a  binomial  coefficient. 


r'  rt(n~r)l‘ 

an^i(q)  (i  "  0,1,*.* n-1)  i  ■»  the  coefficienta  of  the  polynomial  defined  by 


n-1 


l  «n  < 

i-0  n,i 


— - — -  n  (x»q).  v  t  ” 

Y  (c-1)  nit" 

n 


an"*i  *  0,1,. ..n-1,  j  -  0,...i(n-1-i))t  -  the  coefficient*  of  polynomial  defined  by 


an,i(q)  - .  *<n-i,<n-i-i,/2  UT“  <j 


1,+- 


Given  a  biinfinite  geometric  knot  sequence  t  i  »  (q  )_ct  for  some  q  e  (0,“)  with 


t  *  -  lim  t  ,  I  »  -  (t_,t  ) 

i**" 


we  denote  by 


Nn,i  '  “  ’Vn'V  Cti'ti+l'*”ti+n3  (*‘x,+ 


n-1 


the  corresponding  B- splines  normalized  so  that 
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and  by  S  i  ■  span{N  },  tha  spaca  of  aplinea  of  degree  n-1  with  knots  t. 

lit  1 

We  can  consider  the  projectors  P  t  C(I)  ♦  S_.  _  defined  by  the  condition  that 

R#t 


Pk,rf  "  |  “i  <f)N2k-r,i 
\  (Nr.i'N2k-r,j»*j(f)  ' 


(N 


r,i 


,f ) 


b 

with  (f,g)  s  »  /  f(x)g(x)dx. 

a 

Then  Pk  0  is  the  interpolation  projector  and  Pk  k  the  usual  L2 -projector  onto  Sk  t 
This  paper  is  a  continuation  of  (11*  In  [1]  the  uniform  boundedness  of 

IG 

for  q  e  (0,»)  with  r  «  k,  k-1  was  proved.  Here,  g”1  is  the  inverse  of  corresponding 
B-spline  Gram  matrix.  In  this  paper  we  obtain  the  sharp  lower  bound  for  IG  1 I  .  I.e.  we 
prove  that  for  any  q  t  (0,~)  and  r  -  k,  k-1,  the  inequality 


n2K-1(qr,q> 

n2k-1('qr,q)l 


n2k..(qr»q>  I 

-  >  2k 

n2k-1(-qr,q,i 


1 


holds. 

In  order  to  prove  this,  we  need  some  properties  of  H  (1  »q)  which  were  studied  in  [1] 


and  [2].  Por  the  reader's  convenience  we  copy  some  of  them  as  follows 


I 


Proposition  1  [2]  IMXiq)  satisfies  a  "difference- delay"  aquation 
nQ(X;q)  «  ■  1 

nn+1<X»q)  -  ((1-X)qBnn(q"1X»q)  -  (qn+1-X)iyX»q) ),  n  -  0,1,...  . 


Proposition  2  [1]  The  polyncnial  IIn(Xiq)  satisfies 


_  , ,  ,  ,n-1  -n(n-1)/2_  ,  n, -1  . 

n  (X;q)  -  X  q  (q  X  »q) 

n  n 


(i.D 


The  coefficients  an,±(<X)  can  be  computed  recursively  by 


where 


an+i,i(q)  "  «rif1<<qn+1-q,*‘1>*n<1<<i>  ♦  (q^'^-Ha^j^iq))  ,  (1.2) 


an,0<‘*)  *  -  *n,-l‘q>  “  «n,nlq>  ‘  "  0 


Proposition  3  [1]  The  coefficients  a^^q)  satisfy 


S,lW  - 


{1.3) 


and  for  n  >  2  the  integer  coefficients  a^J  are  symmetric 

n,i 


In  particular 


(1.5) 

*  "•«> 

2.  The  sharp  lower  bound  for  IG^1 1^ 

Before  proving  the  theorem  we  need  to  do  some  preparation. 

Lena  2. 1  The  following  equalities  hold 


*n,i 


l  i  J 


‘(1! 

n,i 


<n-2)  ("-1)  -  | 


*2k-1<q-1) 


k  .  3 

1  ,  .k-1  2 

2k^r  ’  <->  q 


n^.jl-q  iq)  H  jk(W)  (k-1) 


j-0  <°2k“1»  J  ’  B2k-1.))q:I 


with 


(j)  .  r  j-i(2i-D)  ( j-i(2i+1 ) )  \ 

2k-1,j*  ’  2k-1,k-1  i"1  (  2k-1,k-1-2i  2k-1,k-1-2i' 

(j)  +  7  f  <j-K2i-1»  (j-i(2i*1))1 

B2k-1,j  2k-1,k-2  1£1  *•  2k-1,k-2i  2k-1,k-2-2iJ 


(2.1) 


For  convenience,  here  and  below  we  use 


a(r> 

n,i 


0,  if  r  <  0  or  r  >  i(n-l-i)  as  well  as  i  <  0 


(2.2) 


In  particular 


-4- 


Similarly 


with 


Proof » 


with 


n2k-2<-qlt~1,q) 


2k -2 


.  ,  i  <k-1)<3k-4)  <k-1)(k-2) 

, '“a-a.j  ‘  e2k-a,j> 


(-)  q' 


«  Jj>  *  7  f.tj-Kai-D)  .  .(j-l(2i-M)) 

a2k-2,j  °2k-2,  j  1  2k-2,k-2  2k-2,k-1-2i  2k-2,k-2-21 

-  (j-jid+D) 
i-0  *2)t“2'lc"2"1 


By  the  definition  of  an>  ^q)  and  a^j 


n2k-i(-qk'q)  2r2  .  k.  k.i 

— 45 - ^—7  -  l  atqii-q  ) 

Y2k-1(q-1,2k1  1-0 


2k-2  i(2k-2-i) 

-  I  l 

i-0  j-0 


«4+j 


f-i4»  1  a(J) 

2k-1,i 


®i  *  “  2  (2k-1_1><2,c_2"1>  +  ik 


|  k(k-1 ) 

aln  ♦.  -  q 

0<i<2k-2  1 


*1  *  “  *i  k(k_1)  “  \  (k-1-i><k-2-i> 


th«n 


n2k-1  t"ql',q)  f  k(Jc-1)  2k-2  l(2k-2-i) 
. . . q  * 


Y2k-1(q-,) 


2k- 1 


I  I  *-»4  -21,  4q  1 

1-0  j-0  2k  1,1 


*,+d 


f  k(k-1) 


V  (k-^2-i2  .  .k-1-1  7  (J, 


I  I  <->* 

1-0  j-0 


2k-1,k-1-l 


k-1  ‘k-1>2-i2  k.1+ij‘(‘+')+3  (j) 


♦  I  I  (->  q 

1-1  j-0 


2k- 1 ,k+l-1 


,  ,k-1  2  k(k-1» 
(-)  q 


(k-1)2  (j)  k-1  (k-1)1-!2 


I  “2k-1,k-1q3+  I  I 

j-0  '  !»,  jw 

<«  +/1U+1,+jaU>  ) 

2k-1  ,k-1-i  q  2k-1,k-1-i'' 


,  .1,  2  1(1“1)+3 
(-)  (q  t 


v  .  |  k  ( k- 1 )  <k-1)2 
k-1  2  r 

<-)  q  i 

j-0 


(a<i>  -  ,<J>  w  V  ,.<j-i<2i-1)> 

'  2k-l,k-1  ®2k 


-1,k-2 '  1  2k-1,k-1-21 


( j-l(21+1 ) )  _  r  ,  ( J-X21-1 ) )  (j-i(2i+1)) 

2k-1 ,k-1-2i '  ^  *  2k-1,k-21  2k-1,k-2-2l' 


qJ  . 


and  (2.1)  follows. 

By  (2.1)  and  (1.5)  wa  get 


(0) 

r  2k-2 

'2k-1,0  " 

a2k-1 ,k-1 

'  k-1 

a*0’ 

(2k-2 

2k-1 , 0 

2k- 1  ,k-2 

'k-2 

The  sane  kind  of  argument  proves  (2.4). 


6 


An  analogous  argument  gives 


nac-1<qlC,<it)  |  lc(k-D  U-1> 


Y2k-1(q"1> 


2k-1 


L  ,a2k-1,j  +  B2k-1,j)q 


(2.5) 


j-0 


„<qk“1iq>  7  ot-nu*-*)  <k-i)(k-2> 

2k“2  q2  I  (a 


w*-1’ 


2k-2 


j-0 


2k -2, j  +  82k-2,j)q 


(2.6) 


A  straightforward  calculation  starting  with  the  formula  defining  IMXiq)  leads  to  the 
expreseions 


n2k-1(q  ,q)  ( 2k-1 i 

.  ,.2k-1  *■  k  > 

^2k-1  q""1 


k-1 , 


f  k(k-1)  k-2  , 

q‘  n  (1  +  q+...+q  )  (1  +  q+»**+q  ) 

i-1 


fk(k-l)  (k-1)2 

(  k  ]  q  dat"1'1<3 


(2.7) 


n2k-2(qk  1,q)  f2k-21  2 

,2k- 2  '■k-l  J  4 


Y2k-2<q-1> 


i  (k-1 )(3k-4)  k-2  .  , 

H  (1  +  q*-'**+q  ) 
i-1 


„  ,  |  (k-1 ) ( 3k-4)  (k-1) (k-2) 

W  **  lQ  ^2k-2 , jq 


j 


(2.8) 


From  (2.5),  (2.6),  (2.7),  (2.8)  it  is  easy  to  find  the  following  relations 


(T1) 


2k-1,i 


a2k-1,i  *  S2k-1,i 


(2.9) 


-7- 


1 


(2lc_2l  d  -  a  +  B 

'•k-l  ‘  2k-2,i  2k-2,i  2k-2,i 


d2k-1,i  “  ^  d2k-2,i-j 


k-2 

d  =  T  d 

2k-2,i  2k-3,i-j 


d2k-1  i  *"  d  2 

2k-1,(k-1)  -i 


d2k-2,i  ”  d2k-2,(k-1)(k-2)-i 


d  .  i  «  a  :  =  0  if  i<0 
n,i  n,i 


Lemma  2 . 2  The  following  equality  holds 


»  -  I  »  •  +  I 

n'1  i-t-i  n‘1'1  r»i+i-n+1 


Proof  By  (1.2)  and  the  definitin  of  a 


,  .  (n-i) (n-1-i)/2  1(lr1“1)  (*>  t 

an,l<q)  q  Vi  q 


(2.  10) 


(2.11) 


a  ,(q)  “  q”  1  i<1  +  q*-,**+q1>a  ,  .  iq)  +  (1  +  q+*«.+qn  i_1)  a 

n / i  n- 1,1 


(n-i)(n-1-i)/2 


i(n-l-i)  j  ,  ,  1 

i  t  l  *ir ,  *  i 

j-0  r-j-i  n  r-J+i-n+1  1,1  1 


n-1 , i-1 
(r> 


(q) 


(2. 16) 


Comparing  with  (*),  !2.15)  follows. 


Corollary  2. 1  The  following  inequalities 

4  «  [j  i(n-i-l)],  0  <  L  <  n-1  (2.17) 


0  <  4  <  i  (n-1-i) ,  i  <  (2.  18) 

hold. 

Proof  We  use  mathematical  induction  to  prove  (2.17),  (2.18).  Suppose  for  n-1  (2.17), 

(2.18)  hold.  Using  (2.15), 


n,i 


>  a 


U-1) 

n,i 


for 


and 


a(1>  >  aU>  , 
n,l  n,i-1 


for 


(4)  (i-1) 

,i  ”  n,i 


(a 


(i) 
n-1  ,i 


(i-i-1) 

n-1,i 


)  +  ( 


aU) 
n-1, i-1 


(i+i-n), 
an-1,i-1'  * 


Since  (1.4) 


(4)  (i(n-l-i)-i) 

an,i  3n,i 


and 


4  <  —  i(n-l-i)  as  well  as  0  t  i  <  n-1 


We  get 


-9- 


(1)  U-i-1)  it  ,  t  1  4, 

n-1,1  *  an-1,i  '  lf  1  *  2  1(  2  l) 

U>  .  [2  l<n-2-i)"  21  „  (t-i-1 )  1  „  1  ,, 

a  .  .  ^  ^  •  .  j  »  ^  t  i(n-2-i)  K  t  ^  •  i(n-l-i) 

n  n  i,i  n- i # 1  2  2 


a(A*  >  .“-V>  for  0  <  A  <  ~  i(n-l-i) 
n- 1 1  x  n- 1 ,  x  2 


The  same  kind  of  argument  shows 


(A)  (A+l-n) 

“n-1,1-1  n-1,1-1 


Now,  we  bring  the  induction  hypothesis  to  the  next  level  and  (2,17)  is  proved  since  it 
obviously  holds  for  n  <•  2, 

In  order  to  prove  (2.18),  it  is  enough  to  prove  (2.18),  only  for 
0  <  A  <  j  (i-1)(n-i)  because  of  (2.17)  and  (1.4).  By  (2.15) 


r  /  (r)  (r)  ^  .  r  r  _ ( r )  (r) 


<*>  -  a(t>  .  Y  fa11'  -  a,r'  1  +  Y  fa11'  -  a11'  1 

an,i  an,i-1  rJ;i+1  l  n-l.i  n-1.i-lJ  '-1.1-1  ~1.*-*J 

.  fa(4-l>  -  aU+1_nM 
^  n-1,i  n-1,i-2l  * 


By  induction  hypthesis  and  (2.17),  we  know 


aU)  >  a(r) 
n-1,i-1  n-1,i-2 


(r)  (r)  . 

an-1,i  "  n-1,1-1  0 


as  wall  as 


(t-l)  (*-±)  <i+l-n))  (t+i-n) 

n-1,1  n-1,1-1  n-1,1-1  a- 1,1-2  * 


Therefore  (2.18)  holds  for  n  and  so  (2.18)  is  proved,  since  it  Is  obviously  right  for 


Lemma  2.3  The  following  equalities 


°2k-1 , j  “  2 


3  -  j-i(2i-1)  .  . 

1  a  r  +  l  {  Y  a(r) 

jHfl  2k-2'k-1  i*1  r"  j-k+1-i(2i-3)  2k-2,)c-1-2i 


+  (r)  +  (r, 

r-j-k+1-i(2i+1)  2k-2'k“2"21  r-j-k+1-i(2i-1 )  2k-2'k'1_2i 


j-i(21+1) 

+  .  .  J _  a2k-2,k-2-2i^  (2‘ 


r-3-k+1-l(21+3) 


J  .  .  )  a<r)  +  7  ,  (r, 

i-j-k+1  2k-2,i  r-j-k+1  2k’2>k"2  i_!  r- j-k+1-i(21-1)  2k-2'k-1' 


(2.20) 


3-1(21+1 ) 

+  Y  Jr>  1 

L  2k-2 . k-2-2 i '  ' 

r-3-k+1-i(21+1)  ' 


Proof  From  (2.1),  (2.4)  and  Lemma  2.2,  by  straightforward  calculations,  (2.19)  and  (2.20) 


can  be  proved. 


Lemma  2.4  For  3  <  (j  (k-1)2),  the  inequality 


a  <  a  *  d 

2k- 1 ,  j  2k- 1,0  2k- 1 ,  j 


Proof  Sine* 


til  *  d 

2k-1,0  2k- 


3  "  (^)  1  d2k-2,t 

*3  K  1  t-j-k+1  ' 


“  2  I 


i-j-k+1 


2k-2,i 


In  order  to  prove  Lemma  2.4  it  ia  enough  to  show 


4  2  t-jLi aa-2'1 


(2.21) 


(2.22) 


From  (2.19)  and  (2.20) 


J 

i«A+i  a2k-2'4  ‘  a2k“1'i 

-  j-k-i(21-3) 

"  I  (  l  JT) 


j-l(2i-1> 

2V-2  .k-1-2i  +  I 

i-1  r-j-k+1-l(2i-1)  r-j+1-i{2i*1) 


s(r) 

2k-2,k-1-2i 


3-i(2i-1,  _(l) 

r-j+1-l(2i+1) 


j-k-i(21+1)  (r>  , 

2k-2,k-2-2i  2k-2,k-2-2i-1 


r-3-k+1-i(2i+3) 

-  j-i(2i-1) 

y  f  y  fa(r)  -  a<r>  1 

1-1  r— j+1-i(2i+1 )  1  2k-2,k-1-2i  2k-2,K-2-2iJ 


J-k-i(2i-3) 


+  y  r.w  -  ,tr"41)  n 

r-JHc+1-l(2i+1)  1  2k-2.k-1.2i  2k-2 .k-2-21  * 


Because  of  (2.17),  (2.18)  and  (1.4),  Lemma  2.4  ia  verified. 

■ 

After  the  preceding  preparations  now  it  is  time  to  prove  the  following  theorem. 
Theorem  For  any  q  e  (0,«)  and  r  -  k-1,k,  the  inequality 


IG 


-1. 


n2k-1(q  ,q) 


n2k-l<V'q) 


>  2k- 1 


(3.0) 
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hold*/  and  lim  IG_1I  ■  2k- 1 . 
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Proof  Bacauaa  of  the  symmetry  of  n  (X»q),  w«  can  raatrict  our  dlacussin  to  the  caaa 

n 

q  e 

Since  (2.1),  (2.3),  (2.7),  Lemma  2.4  and  (1.1),  we  know 


n2k-1(qt"1,q) 

|n2k-1(qlt’q)  1 

_  k-1 

k  i 

n2K-1(-q  ,q) 

*n2k-1<“q  ,q)* 

,v  ,  (k-D 

(T1)  I 


i-0 


2k-1,i^ 


Ot-1)  ...  . 

£  (2a2k-1,i  "  (  k  ^  d2k-1,i)q 


(k-1) 


r  2k-1  s  .  i 

*  k  J  J  d2k-1,lq 
1<*0 


(k-1 ) 


?  (2a  •  d  -  f2k-1]  d 

i“0  '  2k-1,0  a2k-1,i  1  k  >  2k-1,J 


-  2k -1  , 


and  refer  to  [1]  for  equality. 

the  proof  of  the  theorem  reliea  on  Lenma.  2.4  mainly.  In  order  to  prove  the 
monotonicity  of 


n2k-1(qlC,q)  | 
n2k-1<_qlt'q)  ^ 


for  g  e  (0,»)  , 


it  is  sufficient  to  prove  the  stronger  inequality 


3 

j 

< 

I  ! 

t  i 

w#  fail  to  prove  this  inequality .  But  numerical  result*  (see  appendix)  show  the  inequality 
is  true  at  least  for  n  <  9. 
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